Search for: All records

ABSTRACT The origin of the ‘seeds’ of supermassive black holes (BHs) continues to be a puzzle, as it is currently unclear if the imprints of early seed formation could survive to today. We examine the signatures of seeding in the local Universe using five $$[18~\mathrm{Mpc}]^3$$BRAHMA simulation boxes run to $z=0$. They initialize $$1.5\times 10^5~\rm {M}_{\odot }$$ BHs using different seeding models. The first four boxes initialize BHs as heavy seeds using criteria that depend on dense and metal-poor gas, Lyman–Werner radiation, gas spin, and environmental richness. The fifth box initializes BHs as descendants of lower mass seeds ($$\sim 10^3~\rm {M}_{\odot }$$) using a new stochastic seed model built in our previous work. In our simulations, we find that the abundances and properties of $$\sim 10^5-10^6~\rm {M}_{\odot }$$ local BHs hosted in $$M_*\lesssim 10^{9}~\rm {M}_{\odot }$$ dwarf galaxies, are sensitive to the assumed seeding criteria. This is for two reasons: (1) there is a substantial population of local $$\sim 10^5~\rm {M}_{\odot }$$ BHs that are ungrown relics of early seeds from $$z\sim 5-10$$; (2) BH growth up to $$\sim 10^6~\rm {M}_{\odot }$$ is dominated by mergers in our simulations all the way down to $$z\sim 0$$. As the contribution from gas accretion increases, the signatures of seeding start to weaken in more massive $$\gtrsim 10^6~\rm {M}_{\odot }$$ BHs, and they are erased for $$\gtrsim 10^7~\rm {M}_{\odot }$$ BHs. The different seed models explored here predict abundances of local $$\sim 10^6~\rm {M}_{\odot }$$ BHs ranging from $$\sim 0.01-0.05~\mathrm{Mpc}^{-3}$$ with occupation fractions of $$\sim 20-100~{{\ \rm per\ cent}}$$ for $$M_*\sim 10^{9}~\rm {M}_{\odot }$$ galaxies. These results highlight the potential for placing constraints on seeding models using local $$\sim 10^5-10^6~\rm {M}_{\odot }$$ BHs hosted in dwarf galaxies. Since merger dynamics and accretion physics impact the persistence of seeding signatures, and both high and low mass seed models can produce similar local BH populations, disentangling their roles will require combining high and low redshift constraints. 

Abstract Detecting the first generation of stars, Population III (Pop III), has been a long-standing goal in astrophysics, yet they remain elusive even in the JWST era. Here we present a novel NIRCam-based selection method for Pop III galaxies, and carefully validate it through completeness and contamination simulations. We systematically search ≃ 500 arcmin²across JWST legacy fields for Pop III candidates, including GLIMPSE, which, assisted by gravitational lensing, has produced JWST’s deepest NIRCam imaging thus far. We discover one promising Pop III galaxy candidate (GLIMPSE-16043) at $z=6.50_{-0.24}^{+0.03}$, a moderately lensed galaxy ( $\mu =2.9_{-0.2}^{+0.1}$) with an intrinsic UV magnitude of $M_{\mathrm{UV}}=-15.89_{-0.14}^{+0.12}$. It exhibits key Pop III features: strong Hαemission (rest-frame EW 2810 ± 550 Å); a Balmer jump; no dust (UV slopeβ = −2.34 ± 0.36); and undetectable metal lines (e.g., [Oiii]; [Oiii]/Hβ < 0.44), implying a gas-phase metallicity ofZ_gas/Z_⊙ < 0.5%. These properties indicate the presence of a nascent, metal-deficient young stellar population (<5 Myr) with a stellar mass of ≃10⁵M_⊙. Intriguingly, this source deviates significantly from the extrapolated UV–metallicity relation derived from recent JWST observations atz= 4–10, consistent with UV enhancement by a top-heavy Pop III initial mass function or the presence of an extremely metal-poor active galactic nucleus. We also derive the first observational constraints on the Pop III UV luminosity function atz ≃ 6–7. The volume density of GLIMPSE-16043 (≈10⁻⁴cMpc⁻³) is in excellent agreement with theoretical predictions, independently reinforcing its plausibility. This study demonstrates the power of our novel NIRCam method to finally reveal distant galaxies even more pristine than the Milky Way’s most metal-poor satellites, thereby promising to bring us closer to the first generation of stars than we have ever been before. 

Abstract We study the black hole mass–host galaxy stellar mass relation,M_BH–M_*, of a sample ofz< 4 optically variable active galactic nuclei (AGNs) in the COSMOS field. The parent sample of 491 COSMOS AGNs were identified by optical variability from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) program. Using publicly available catalogs and spectra, we consolidate their spectroscopic redshifts and estimate virial black hole masses using broad-line widths and luminosities. We show that variability searches with deep, high-precision photometry like the HSC-SSP can identity AGNs in low-mass galaxies up toz∼ 1. However, their black holes are more massive given their host galaxy stellar masses than predicted by the local relation for active galaxies. We report thatz∼ 0.5–4 variability-selected AGNs are meanwhile more consistent with theM_BH–M_*relation for local inactive early-type galaxies. This result is in agreement with most previous studies of theM_BH–M_*relation at similar redshifts and indicates that AGNs selected from variability are not intrinsically different from the broad-line Type 1 AGN population at similar luminosities. Our results demonstrate the need for robust black hole and stellar mass estimates for intermediate-mass black hole candidates in low-mass galaxies at similar redshifts to anchor this scaling relation. Assuming that these results do not reflect a selection bias, they appear to be consistent with self-regulated feedback models wherein the central black hole and stars in galaxies grow in tandem. 

Abstract We make an in-depth analysis of different active galactic nuclei (AGN) jet models’ signatures, inducing quiescence in galaxies with a halo mass of 10¹²M_⊙. Three jet models, including cosmic-ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. Each of our simulations is idealized isolated galaxy simulations with AGN jet powers that are constant in time and generated using GIZMO and with FIRE stellar feedback. We examine the distribution of Mgii, Ovi, and Oviiiions, alongside gas temperature and density profiles. Low-energy ions, like Mgii, concentrate in the interstellar medium (ISM), while higher energy ions, e.g., Oviii, prevail at the AGN jet cocoon’s edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic-ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic-ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in Oviand Oviiiis demonstrated in the ISM-to-circumgalactic medium (CGM) transition (between 10 and 30 kpc), growing smoothly toward the CGM at greater distances. 

Abstract Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofR_B≈ 2 × 10⁵GM_•/c², whereM_•is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r⁻¹rather thanr^−3/2and the accretion rate $\dot{M}$is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{\mathrm{B}}$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^2\sim 5\times {10}^{-5}{\dot{M}}_{\mathrm{B}}{c}^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback. 

Abstract The recent Chandra-JWST discovery of a quasar in thez≈ 10.1 galaxy UHZ1 reveals that accreting supermassive black holes were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity ofL_bol∼ 5 × 10⁴⁵erg s⁻¹, which corresponds to an estimated black hole (BH) mass of ∼4 × 10⁷M_⊙, assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, overmassive black hole galaxies (OBGs) by Natarajan et al., that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Given the excellent agreement between the observed multiwavelength properties of UHZ1 and theoretical model template predictions, we suggest that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux, which is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift ofz≈ 10.1, as predicted for the transient OBG stage (9 <z< 12); the amplitude and shape of the detected JWST spectral energy distribution (SED) between 1 and 5μm, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1, which is suggestive of a recent merge and is also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe. 

While supermassive black-hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event-horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black-hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes. 

Abstract Supermassive black hole binary systems (SMBHBs) should be the most powerful sources of gravitational waves (GWs) in the universe. Once pulsar timing arrays (PTAs) detect the stochastic GW background from their cosmic merger history, searching for individually resolvable binaries will take on new importance. Since these individual SMBHBs are expected to be rare, here we explore how strong gravitational lensing can act as a tool for increasing their detection prospects by magnifying fainter sources and bringing them into view. Unlike for electromagnetic waves, when the geometric optics limit is nearly always valid, for GWs the wave-diffraction-interference effects can become important when the wavelength of the GWs is larger than the Schwarzchild radius of the lens, i.e., $M_{\mathrm{lens}}\sim {10}^8{\left(\frac{f}{\mathrm{mHz}}\right)}^{-1}{M}_{\odot }$. For the GW frequency range explored in this work, the geometric optics limit holds. We investigate GW signals from SMBHBs that might be detectable with current and future PTAs under the assumption that quasars serve as bright beacons that signal a recent merger. Using the black hole mass function derived from quasars and a physically motivated magnification distribution, we expect to detect a few strongly lensed binary systems out toz≈ 2. Additionally, for a range of fixed magnifications 2 ≤μ≤ 100, strong lensing adds up to ∼30 more detectable binaries for PTAs. Finally, we investigate the possibility of observing both time-delayed electromagnetic signals and GW signals from these strongly lensed binary systems—that will provide us with unprecedented multi-messenger insights into their orbital evolution. 

Abstract Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ∼10–100 Hz, over a characteristic jet activity timescale ∼10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ∼10 events yr −1 . The detection rate drops to ∼1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events. 

We present estimates for the number of supermassive black holes (SMBHs) for which the next-generation Event Horizon Telescope (ngEHT) can identify the black hole “shadow”, along with estimates for how many black hole masses and spins the ngEHT can expect to constrain using measurements of horizon-resolved emission structure. Building on prior theoretical studies of SMBH accretion flows and analyses carried out by the Event Horizon Telescope (EHT) collaboration, we construct a simple geometric model for the polarized emission structure around a black hole, and we associate parameters of this model with the three physical quantities of interest. We generate a large number of realistic synthetic ngEHT datasets across different assumed source sizes and flux densities, and we estimate the precision with which our defined proxies for physical parameters could be measured from these datasets. Under April weather conditions and using an observing frequency of 230 GHz, we predict that a “Phase 1” ngEHT can potentially measure ∼50 black hole masses, ∼30 black hole spins, and ∼7 black hole shadows across the entire sky.